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Abstract. We present an extended analysis of the relation between radio surface brightness and diameter - the 
so-called E — D relation for planetary nebulae (PNe). We revise our previous derivation of the theoretical E — D 
relation for the evolution of bremsstrahlung surface brightness in order to include the influence of the fast wind 
from the central star. Different theoretical forms are derived: E oc D^ 1 for the first and second phases of evolution 
and E oc D~ 3 for the final stage of evolution. Also, we analyzed several different Galactic PN samples. All samples 
are influenced by severe selection effects, but Malmquist bias seems to be less influential here than in the supernova 
remnant (SNR) samples. We derived empirical E — D relations for 27 sample sets using 6 updated PN papers 
from which an additional 21 new sets were extracted. Twenty four of these have a trivial form of f3 m 2. However, 
we obtain one empirical E — D relation that may be useful for determining distances to PNe. This relation is 
obtained by extracting a recent nearby (< 1 kpc) Galactic PN sample. 

Key words, planetary nebulae: general - radio continuum: ISM - methods: analytical - methods: statistical - 
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1. Introduction 

Planetary nebulae (PNe) are usually identified by their op- 
tical spectra, consisting mainly of recombination lines of 
hydrogen, helium and collisionally excited light elements. 
They are composed of gaseous shells ionized by a hot cen- 
■ tral star, which allow for the necessary conditions to create 
their spectra. PNe were first discovered over 200 years ago. 
Today, there are more than 2500 confirmed Galactic PNe 
with total number estimates at least ten times higher. PNe 
emission from radio to X-ray have been detected over the 
past 30 years. Several hundreds of these PNe have been 
observed in radio-continuum alone. Their number far out- 
rank the known Galactic supernova remnant (SNR) pop- 
ulation (approximately 250). 

The relation between radio surface brightness and di- 
ameter for SNRs, the so-called S — D relation, has been the 
subject of extensive discussions over the past forty years. 
This £ — D relation has a power law form: 



Observational improvements, including the increased 
use of radio interferometers, account for the high number 
of resolved radio PNe over the past two decades. The E— D 
relation for PNe has only been sporadically discussed (for 



example, lAmnuel et al. 1984)) until Urosevic et al. (2007) 
(hereafter Paper I). 

The statistical analysis of the Galactic PN 
distances 



was presented in Daub (1982) 



Cahn, Kaler & Stanghellini (1992) Using 



and in 
radio data, 

the primary statistical method for determining distance 
was based on the correlation between PN radius and 
brightness temperature, known as the R — Tb relation^. 

derived 



Van de Steene fc Zijlstra (1995) 



an empirical 
felt to be located 



R — Tb relation for 131 Galactic PNe 
near the Galactic center, at approximately the same 
distance. Additionally, a smaller sample of 23 (non- 
Galactic bulge) PNe were studied having well-determined 
distances. The same calibration relations were obtained in 



(1) 



both cases. Zhang (1995) studied the statistical distance 



where £ is the radio surface brightness, D is diameter and 
P the slope of the log-log plot. 



scale for PNe based on another sample of 132 PNe having 
well-determined individual distances. He found that 
the derived distance scale was in good agreement with 
those derived by Van de Steene & Zijlstra (1995) More 
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This relation is essentially equivalent to the D — E relation. 
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recently, |Phillips (2002)] (hereafter Ph02) derived R-T h 
relations for 44 nearby PNe whose distances are less than 
0.7 kpc. Thus, different samples of Galactic PNe with 
known distances were defined in the three previously cited 
papers. All of these derived R — T h relations were applied 
to PNe with unknown distances. We also note that 



a recent paper by Stanghellini, Shaw & Villaver (2008) 



(hereafter SSV) gives an updated reliable sample of 
individual Galactic PN distances. 

The statistical analysis of Galactic radio PNe is less 
developed than those of Galactic Therefore, selec- 

tion effects are felt to greatly influence Galactic radio PN 
samples causing the statistical determination of distance 
to be highly uncertain. 

In the present paper, we use the methodology for the 
statistical analysis of SNRs to improve that for Galactic 
radio PNe. Primary objectives include the following: 

(i) Development of an improved form of the theoretical 
£ — D relation for PNe (compared to the one derived 
in Paper I) by incorporating the interaction between 
asymptotic-giant-branch (AGB) star wind and the 
fast wind from the central star. 

(ii) Derivation of improved empirical £ — D relations 
using current and previous samples to determine if 
these are all under the influence of selection effects. 

(iii) Validation of £ — D relations for the determination 
of radio PN distance. 



radiation mechanism used for derivation of the theoretical 
£ — D relation is thermal bremsstrahlung, which is the 
basic process of radio-continuum emission in H II regions. 
We begin our derivation using the following notation: 

R s - shell radius, 

R - outer radius of the shell, 

Ri - inner radius of the shell, 

V s - velocity of the shell, 

V - AGB wind velocity, 

v - fast wind velocity, 

M - mass loss in the AGB phase, 

m - fast wind mass loss, 

M s - mass in the shell, 

p s - density of the shell, 

t = - the moment when the AGB wind stops, 
t = t - the moment when the fast wind starts. 



2. Dynamics of PNe 

In considering the dynamics of PNe, we adopt the in- 



teracting stellar winds model (ISW; | Zhang fc Kwok 1993 



IKwok 1994[ and references therein). This model implies 
that a nebular shell results from the sweeping up of cir- 
cumstellar material moving within the AGB star wind by 
a fast wind emanating from the PN central star. The speed 
of the fast wind is 100 times that of the AGB wind and acts 
as a " snowplow" , piling up matter into a high density shell. 
This interaction produces two shocks. The inner shock is 
formed at the interaction between central star wind and 
shell. The outer one is situated in the AGB wind just out- 
side the shell. Most of the volume interior of the shell is 
made up of shocked central-star wind at a high tempera- 
ture (millions of degrees). The outer shock is likely to be 
isothermal - similar to the shocks in the radiative phase 
of SNR evolution, where a dense shell is formed. 

The dependence between gas density and radius of a 
PN is necessary for deriving the theoretical £ — D rela- 
tion for PNe. This dependence was shown in Paper I for a 
model consisting of one steady wind emanating from the 
central star. A trend of decreasing radio surface bright- 
ness with increasing diameter was shown in Paper I. The 



2 The radio samples of Galactic PNe are richer in terms of 
the number of objects, but incompleteness is drastically higher 
than in SNR samples - less than 5% of total number of PNe 
can be detected by radio surveys (|Kwok 1994|) . 



2.1. The momentum conserving phase 

If we assume momentum conservation during the interac- 
tion, then: 



P, = M.V. = P 



s v s — - 1 AGB ~r A FW 1 



Ps — Page ~\~ Pfw 

4nr 2 p AGB Vdr + / 4irr 2 p FW vdr 

Vt JRi 



(2) 



4nr 2 



M 



■Vdr 



v(t-r) 



m 

4irr - — 75— vdr 



ATtr 2 V J R . 4Trr' 2 v 

= M(R -Vt)+rh(v(t-T)-Ri), 

while the mass in the shell is: 



(3) 



..2 . 



M s = M AG _ 
Ro 

^■nr z p AGB dr 

vt 

R ° 2 M , 
inr - — — dr 

vt 4:TTr 2 V 



ti(t-r) 



4-Kr 2 p FW dr 



Ri 



ti(t-r) 



47rr" 



Airr 2 v 



dr 



= ^-(R -Vt) + ™(v(t~T)~R t ). 
V v 



(4) 
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If, additionally, the shell is thin: i?,; — > R a = R s then we 
may write: 



P s = (M -m)R s - MVt + mv(t - t), 
,M m N 



M, 



V v 



)R S -Mt + m{t-r) 



(5) 
(G) 



(see Equation 6 in|Zhang fc Kwok 1993[ ). 

If V s = const., the shell is expanding following the law: 



_ VVT , VT 

R, = - + V 3 (t- 



v-V 



v-V 



(7) 



which may be reduced to Equation (5) in 
Zhang fc Kwok (1993)|if v > V: 



R s ^Vr + V s (t-T). 
If t > r then, 
R s » V.t, 

and the shell velocity can be expressed as: 



p s ( M - m) R s + {rhv - MV) % 



Ms (m_^ Rs + {i}i _m)^ 
Solving this equation we obtain: 



(8) 



(9) 



(10) 



V s 



M m 
V v 



(11) 



which is Equation (4) in Zhang & Kwok (1993) 



The average density inside the shell can be expressed 



Ps ~ ^.fRl ' 



(12) 



where / is the volume filling factor. From Equations (6) 
and (9) we see that if t ^> r: 



Ms 



,M rhs m — M 



Vs 



R s 



If M > m, and v > V, then M s ra {f + ^ 
M 



(13) 
^)R S and 



(14) 



Essentially, we have perturbed the AGB wind 
(Pagb = v ). From Equation (11) we es- 

AGB 

timate Vg ra 15 km s _1 using typical values of 
M = 1Q- 5 Mqyt- 1 , ra = 10" 8 Afoyr" 1 , V = 
10 km s _1 and v = 2000 km s~ . We obtain the depen- 
dancc p s oc R^ 1 using Equations (12) and (13) assuming 
that the shell thickness A = R a — i?; = const, and 
/ = 1 - (ft) 3 = ^) 3 « 3A/i? s . 



2.2. The energy conserving phase 

The momentum-conserving phase is perhaps, a reasonable 
description of early PN evolution. However, if a part of the 
energy of the fast wind is transformed into thermal energy, 
then the pressure of the hot bubble (P) will provide ad- 
ditional acceleration to the nebular shell. In this "energy- 
conserving" phase, the equation of motion and the energy 
equation is expressed as: 



M s ^=^RlP~^{Vs-V)\ 



1 4tt , 
P ■ — R 3 

dtV 7 -l 3 8 



1 ■ 2 

—mv 
2 



i-KRtPVs 



(15) 



(16) 



(|Kwok 2000|) . 

The similarity solution method can be used for defining 
the change of PN radius with time. The power-law form 
R s oc t a , is expected for a PN in the energy conserving 
phase of evolution. We can define dimensionless variable 
in form: 



£ = R s t l p™^T>E n 



(17) 



where E is the power of fast wind which injects energy 
(through the inner shock) into the hot bubble. After sim- 
ilarity analysis, the following expression is obtained: 



Rs=^rti. 

Pagb 



(18) 



We obtain a R s oc t (a = 1) dependence assuming 
Pagb oc R~ 2 and using Equation (18). Finally, a station- 
ary solution V s = R s = const, can be used for further 
derivations. 

The mass of the compressed AGB wind in the shell, 
using V^const., is: 



Ms = M(±-±-h 



(19) 



If for /, we assume A = const., then Equation (12) implies 
p s oc R- 1 . 

2.3. The final phase 

We only have AGB wind momentum and mass conserva- 
tion in the final phase of evolution when the pressure of 
the hot bubble, P ->■ 0: 



(20) 



P s = / 4xr 2 p AGB Vdr + P' 

JR' 

= M(Rs-R') +P' 



M s = / 4nr 2 p AGB dr + M' 

JR 1 

M 



V 



(Rs - R') + M', 



(21) 
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where P' and M' are the momentum and mass acquired 
by the time the shell radius reaches some value R' . For 
the shell velocity we now have: 



M(R S -R')+P' 



M R, - R' 



M'V 



V 



v, 



(22) 



as R s — > oo. 

During the early phases, the PN shell tends to the con- 
tact discontinuity, while in the last phase, it will be just 
behind the isothermal forward shock. If we apply isother- 
mal shock jump conditions, we can estimate: 



p a = M p AGB , 



(23) 



where M = V s /c s is the Mach number and c s , the isother- 
mal sound speed. The density jump is not limited to a 
maximum value of 4, as in an adiabatic shock. If V s ~ 
const., then p oc Rj 2 behind the shock. 

3. Theoretical £ — D relation for PNe 

Assuming a thermal brcmsstrahlung mechanism is respon- 
sible for radiation of H n regions at radio wavelengths, we 
may write the volume cmissivity e v of a PN as: 



£j,[ergs s" 1 cm 3 Hz l ] oc n 2 T 1//2 , 



(24) 



(|Rohlfs fe Wilson 1996|) where n is the volume density and 
T the electron temperature. 

The surface brightness can be expressed as: 



(25) 



where L v is luminosity and D = 2R S , the diameter of the 
PN. We combine Equations (24) and (25), use the derived 
dependence from Section 2 such that n oc p s oc D , and 
assume that temperature is constant [f| This radio surface 
brightness to diameter relation has a different form for the 
first and second phase of evolution (in comparison to the 
earlier result given in Paper I): 



£ oc n 2 T- 1/2 D oc D' 1 . 



(26) 



Generally, we may write n oc D~ x , where x ^ 1 (if A = 
const.). 

We do not expect the temperature to be strictly con- 
stant throughout the nebula. There are temperature gra- 
dients in PNe arising from radiation hardening. More en- 
ergetic photons will travel further and when they are ab- 
sorbed by the PN, they will impart a greater kinetic en- 
ergy to the ions thereby producing a higher temperature. 
This would only slightly increase the value for (3 (+0.1 
approximately), although the T e — D dependence is not 
quite of the power-law form (see numerical model results 
from |Evans fe Dopita 1985| . 

The value (3 = 1 under the above conditions (including 
A = const.) is only a theoretical lower limit; the £ — D 



relation could be steeper due to the these effects, as can be 
seen from Equation (26). The relation given here is derived 
under the assumption of spherical symmetric expansion. 
However, we know that the large number of PNe show a 
high asymmetry and can have a very complex morphol- 
ogy. This may explain why the theoretical £ — D relation 
is not necessary connected to the empirical relations that 
we discuss in Section 4. Even if we had a spherically sym- 
metric PNe with S oc D -1 , the constant of proportionality 
would be different for different phases of evolution, as can 
be inferred from Section 2. Moreover, parameters govern- 
ing the evolution (V, v, M, m) are not the same for all 
PNe, so no unique empirical relation for PNe is feasible. 

In the final stage of evolution, after the fast wind has 
ceased, the isothermal shock wave continuously perturbs 
AGB wind. Nonetheless, the surface brightness steeply de- 
clines, as demonstrated in Paper I, because the mechanical 
energy input (acceleration) from the central star no longer 
exists. Using p AG B oc Rj 2 , the £ — D relation in the final 
stage of PN evolution can be expressed as: 



£ oc D~ 



4. The empirical £ — D relation for PNe 



(27) 



The most important prerequisite for defining a proper em- 
pirical £ — D relation is the extraction of a valid sample 
of PNe. This sample should consist of calibrators having 
well determined distances. 

The distances to the calibrators have to be determined 
by reliable methods. These may include trigonometric or 
spectroscopic parallaxes of PN central stars and meth- 
ods based on nebulae expansion. Although limited by cur- 
rent technology, trigonometric parallax is the most direct 
method of measuring distances. 

All extracted PN samples suffer from severe selection 
effects caused by limitations in survey sensitivity and res- 
olution. For Galactic samples, Malmquist biacl may be 
the most severe effect; this is similar to the situation for 
Galactic SNR samples (jUrosevic et al. 2005| . 

For the literature extracted PN samples in this paper, 
we used flux densities at 5 GHz. This is because most PNe 
are expected to be optically thin at this frequency (see 
Section 5). Radio diameters would be preferable, however 
for most PNe used here, only optical diameters are avail- 
able. 

We extracted 27 PN sample sets as shown in Tables 1 
and 2. While both detail fit parameters for £ — D and 
L-D, Table 1 does so by author (18 PN sets) and Table 2 
by the specific distance determination method used (9 PN 
sets). 

Data for sample Nos. 1 - 4 in Table 1 are taken from 



SSV and Cahn et al. (1992) We used Equation (7) from 



Van de Steene & Zijlstra (1995) to convert the brightness 



temperatures listed in Ph02 to their 5 GHz flux. PN dis- 
tance ranges (kpc) are listed in Table 1 (e.g. "1.0 < d < 



H II regions are approximately isothermal at T ~ 10 4 K. 



4 Intrinsically bright PNe are favored in any flux limited sur- 
vey because they are sampled from a larger spatial volume. 
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2.0" indicates that all PNe in that set have distances be- 
tween 1 and 2 kpc). Three PNe (He 1-5, HDW 6 and 
HaWe 5) indicated by "f" were removed from some sets 
because of inconsistencies with other data, as we will de- 
tail in Section 5. 



Sample set 15 is taken from Zhang (1995) and consists 
of 132 PNe with well-determined individual distances. 
Set 16 contains 23 PNe with well-determined distances, 
while, set 17 consists of 132 Galactic bulge PNe. Both 



samples were taken from Van de Steene & Zijlstra (1995) 



Set 14 consists of 109 Galactic bulge PNe taken from 
Bensby fc Lundstrom (2001)| Set 18 (88 PNe) is actually 
the Van de Steene & Zijlstra (1995) Galactic bulge sam- 
ple with stronger criteria applied to extract Galactic bulge 



PNe as defined in Bensby & Lundstrom (2001) 



The first eight sample sets listed in Table 2 were ob- 
tained from SSV and Ph02 data. Set 9, USNO-PN, con- 
sists of PNe with well-measured trigonometric parallaxes, 
determined by United States Naval Observatory (USNO) 
ground based observations ([Harris et al. 2007j) . Details of 
this USNO-PN sample data are presented in Table 3. 

The trigonometric parallaxes of PNe determined by 
the Hipparcos mission suffer from large measurement er- 
rors (jHarris et al. 2007| and therefore are of limited prac- 
tical use for our study. The average measurement error in 
Hipparcos trigonometric parallaxes for some of our PNe 
is 57 %. Thus, the USNO-PN sample set is more reliable 
and have more accurate distance measurent^). The £ — D 
relation for the USNO-PN sample has the following form: 



7.39t^°10- 



22 



D 



-2.38±0.56 



(28) 



Corresponding Y, v — D and L v — D diagrams are shown 
in Fig. 1. 

Malmquist bias tends to increase the slope of the £ — D 
relation for SNRs (jUrosevic et al. 20 05). Similar implica- 
tions should be valid for PNe as well. Slopes from Table 1 
suggest that some of our sample sets may suffer from a 
Malmquist bias. Specifically, £ — D slopes derived from 
the SSV and SSV+Ph02 are steeper for higher sampling 
distances (i.e. > 2 kpc). 

In spite of the possibility of Malmquist bias, the ma- 
jority of the £ — D slopes listed Tables 1 and 2, are very 
close to the so-called trivial £ — D relation with (3 = 2 
(jArbutina et al. 2004[) . This can be understood mathe- 
matically from: 



£ oc — oc 



S L(D)/d 2 



D 2 /d 2 



oc D- 2 L{D), 



(29) 



where S is the flux density, f2 is the solid angle occupied by 
the source and d is the source distance. The usefulness of 
the £ — D relation can be tested with the correlation coef- 
ficient tl-d of the L — D relation, which should approach 
-1.0 if the data can be explained by a linear relation. 

It is clear from Equation (29) that the £ — D rela- 
tion has a trivial £ oc D~ 2 form in the absence of L — D 



5 An average error of measured trigonometric parallaxes in 
this sample is 17.2 %. 



correlation. The abundance of low tl-d values as shown 
in Tables 1 and 2, indicate that the dependence between 
£ and D is uncertain. Therefore, the validity of £ — D 
relations is diminished even if free from Malmquist bias. 
Aside from a few potential sample sets, including No. 1 
listed in Table 1, most should not be used when calibrat- 
ing the PN £ — D relation. 

Sample No. 1 in Table 1 is the only set with a re- 
spectable correlation coefficient, tl-d = —0.67, yet hav- 
ing a low probability of Malmquist biafl This SSV sample 
consists of 13 relatively nearby Galactic PNe, with reliable 
individually determined distances less than 1 kpc. They 
are different in morphological form (e.g. round, elliptical, 
bipolar core, bipolar, etc.). For this sample, we obtain the 
£ — D relation: 



£ = 7.83t^7 10 -22 jD -2.61±0.21 ; 



(30) 



and a corresponding £, y — D and L v — D diagrams are 
shown in Fig. 2. 

We define the fractional error as: 

dssv — ds 



fe = 



'■SSV 



(31) 



in order to obtain an additional estimate of accuracy. 
Independently derived distances dssv are taken from SSV 
and ds is the distance derived (for each individual PN) 
from the £ — D relation. The maximum and average frac- 
tional errors are fe max = 1.08 and /e = 0.35, respectively. 

According to our criteria and the relatively small aver- 
age fractional error, Equation (30) may be used as a cali- 
bration relation for distance determination. Since the PNe 
in the corresponding sample are nearby, the Malmquist 
bias should not be significant. The L — D relation exists 
and the slope (3 is reasonable. However, this slope is still 
most likely a "mixture" of the evolutionary and selection 
effects. Additionally, the sample number is small (only 13 
PNe) and they do not belong to the same morphological 
type. Therefore, all results obtained using Equation (30) 
should be taken with caution as we discuss below. 

5. Discussion 

In Section 2, we derived the theoretical dependence be- 
tween the gas density in the PN shell and its radius with 
respect to the interacting stellar winds model. Next, we 
showed that the slope of the theoretical £ — D relation is 
changed significantly by including these equations. 

The theoretically derived slope (/? « 1) is shallower 
than slopes from the empirical relations given in Tables 
1 and 2. This discrepancy may possibly be explained by 
poor quality of Galactic PN samples or by incorrect as- 
sumptions used for the derivation of the theoretical rela- 
tion. Due to variations in the power law density distri- 
bution derived here and the approximately constant tem- 
perature of the expanding PN envelope, the theoretical 



6 The correlation coefficients for sets 3 and 9 from Table 1 
are the steepest, but Malmquist bias is also likely. All other 
S — D slopes are approximately trivial. 
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Table 1. The results of the £ — D and L — D fits, where the parameters of the fit (/? and a, respectively) and the 
correlation coefficient r are given for each sample. The number of PNe in the sample is given in the last column (N) . SSV 



stands for Stanghellini, Shaw & Villaver (2008) B&L for Bensby & Lundstrom (2001) and Ph02 for Phillips (2002) 



No. 


Sample" 




T 

' E-D 


a L-D 


r L-D 


N 


n i 


do V a <1.U 


O fil -I-O 01 


n o7 
— u.y * 




— U.o / 




02 


SSV 1 0< d <1 


1 89±0 43 


—0.78 


-0 1 1 +0 

l_/ . _L J 1 KJ ■ 


0.08 


15 


03 


SSV d >2.0 


3.96±0.62 


-0.91 


1.97±0.61 


-0.73 


11 


04 


SSV 


2.56±0.22 


-0.89 


0.56±0.22 


-0.40 


39 


05 


SSV + Ph02 d <0.4 


1.47±0.40 


-0.74 


-0.53±0.40 


0.38 


13 


uu 


oev _i_ pv.no n A <^ d <C\ P, 


o i s+n oo 


n qi 


n i s+n oo 


n 1 7 


O'i 


07 


SSV + Ph02 0.6< d <1.0 


2.37±0.20 


-0.95 


0.37±0.20 


-0.43 


18 


08 


SSV + Ph02 1.0< d <2.0 


2.64±0.40 


-0.88 


0.64±0.40 


-0.41 


15 


00 


SSV + Ph02 d >2.0 


4.44±0.47 


-0.96 


2.44±0.47 


-0.90 


9 


10 


SSV + Ph02 


2.40±0.17 


-0.86 


0.40±0.17 


-0.27 


78 


11 


SSV + Ph02 d <0.4 f 


2.31±0.38 


-0.90 


0.31±0.38 


-0.26 


11 


12 


SSV + Ph02 0.4< d <0.6 t 


2.33±0.15 


-0.96 


0.33±0.15 


-0.44 


22 


13 


SSV + Ph02 f 


2.57±0.14 


-0.92 


0.57±0.14 


-0.45 


75 


14 


B&L 


2.31±0.11 


-0.90 


0.31±0.11 


-0.26 


109 


15 

16 
17 
18 


Zhang (1995) 


2.17±0.14 


-0.80 

-0.84 
-0.92 
-0.94 


0.17±0.14 

0.41±0.34 
0.24±0.09 
0.21±0.09 


-0.10 

-0.26 
-0.24 
-0.25 


132 

23 
132 

88 


Van de Steene & Zijlstra (1995) 


2.41±0.34 
2.24±0.09 

(stronger criteria by B&L) 2.21±0.09 


Van de Steene & Zijlstra (1995) 


Van de Steene & Zijlstra (1995) 



"PNe from Ph02 are included in the sample sets (5 - 13) as additional, if they were not extracted by the SSV in their sample; 
d < 1 denotes that a set of PNe have distances less than 1 kpc, 1 < d <2, between 1 and 2 kpc, etc. 
f - 3 PNe (He 1-5, HDW 6 and HaWe 5) are removed (see text). 

Table 2. The results of the £ — D and L — D fits of PN samples defined using the different methods of distance 
determination. The parameters of the fit (/3 and a, respectively) and the correlation coefficient r are given for each 
sample. The number of PNe in the sample is given in the last column (N). 



No. 


Sample" 




T H-D 


a L-D 




N 


1 


Trigonometric Parallax 


1.79±0.22 


-0.90 


-0.21 ±0.22 


0.23 


18 


2 


Trigonometric Parallax f 


2.02±0.18 


-0.95 


0.02 ±0.18 


-0.02 


17 


3 


Geometrical (trig, parall. + expan.) 


2.26±0.26 


-0.86 


0.26±0.26 


-0.19 


30 


4 


Geometrical (trig, parall. + expan.) f 


2.44±0.22 


-0.91 


0.44±0.22 


-0.36 


29 


5 


Gravitational 


1.42±0.37 


-0.66 


-0.58 ±0.37 


0.33 


22 


6 


Gravitational f 


2.53±0.31 


-0.89 


0.53 ±0.31 


-0.38 


20 


7 


Spectroscopic 


2.57±0.36 


-0.91 


0.57 ±0.36 


-0.43 


13 


8 


Extinction 


1.69±0.38 


-0.80 


-0.31 ±0.38 


0.25 


13 


9 


USNO-PN 


2.38±0.56 


-0.90 


0.38 ±0.56 


-0.32 


6 



" For the first eight samples, PNe are extracted from the SSV and Ph02; for USNO-PN - see Table 3. 
Table 3. The data for the USNO PN sample. 



Name 


Trigonometric Parallax 
[mas] 


SsGHz 

[mJy] 


Diameter 

["] 


NGC 7293 (036.1-57.1) 


4.56 ± 0.49 


1292" 


660 6 


NGC 6853 (060.8-03.6) 


3.81 ±0.47 


1325" 


340" 


NGC 6720 (063.1±13.9) 


1.42 ± 0.55 


384" 


60 c 


A 21 (205.1±14.2) 


1.85 ±0.51 


157 d 


550 e 


A 7 (215.5-30.8) 


1.48 ± 0.42 


305" 


760" 


A 24 (217.1+14.7) 


1.92 ± 0.34 


36" 


415 d 



References: ^Cahn et al. (1992)| lO'Dell (1998)] IGeorge et al. (1974)| ^Hua fc Kwok (1999)| falter et al. (1984) 
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slope may be different than given in Equations (26) and 
(27) (derived assuming the shell thickness A sa const.). 
If the filling factor / is constant, the slope fi could be as 
great as 3. We conclude that whether these simply derived 
theoretical relations have the correct forms or not, empir- 
ical relations are under the influence of biases that make 
their slopes different than they would be if one considers 
the evolution of a single hypothetical PN. An empirical 
relation can be useful for estimating distances only under 
strict conditions as discussed below. 

Van de Steenc & Zijlstra (1995) and Phillips (2004) 
showed that the evolution of PNe may not have a linear 
trend in log- log scales. When the radiation of the central 
star begins to ionize the nebula, the radio brightness will 
increase rapidly. However, in this phase of PN evolution, 
the observed incidence of radio nebulae is predicted to be 
low. Therefore, our linear trend of the £ — D form can be 
justified and the quadratic form of the R — lb dependence 



given by Phillips (2004) is not necessary for the statistical 



analysis given here. In addition, different dependencies can 
not be derived from the thermal bremsstrahlung radiation 
formula (Equation 24). 

The subjects of our analysis are optically thin PNe. 
The reason for the low incidence of observed radio PN dur- 
ing the rising phase, discussed above, is that the majority 
of these PNe arc optically thick at 5 GHz. Then, when the 



central star reaches a temperature of 40 000 K, the num- 
ber of ionizing photons becomes almost independent of 
temperature and radio luminosity no longer increases. As 



discussed by Van de Steene & Zijlstra (1995) this causes 



a decrease in surface brightness with a gradual expansion 
of the nebula, which then becomes optically thin and thus 
detectable. To keep our sample pure, we devised a set of 
criteria which excluded three suspected thick PNe in their 
rising phase of evolution (as marked with "f" in Table 1). 

From Tables 1 and 2, we suggest that 24 slopes 
are most likely trivial X — D relations. Therefore, we 
feel that the Malmquist bias is not so severe as in 
the case of Galactic SNR samples, where the slopes 
are significantly steeper. For example, the obtained 
slope for the sample of 132 bulge PNe collected by 



Van de Steene & Zijlstra (1995) takes a relatively shallow 
value of (3 « 2. This is expected result since all 132 PNe are 
located at the approximately same distance. Malmquist 
bias does not exist in samples which contain objects at 
the same distance. 



Both Schneider & Buckley (1996) and at a later time, 
Bensby & Lundstrom (2001) attempted to improve the 
Bulge sample by removing foreground and background 
PNe. Bensby & Lundstrom (2001) defined stronger cri- 



teria to extract a more pure Bulge sample (Set 14 in 
Table 1). We found a trivial £ — D slope for this sam- 




pie. We also applied this stronger criteria directly to the 
Van de Steene & Zijlstra (1995) sample. The change in 
slope for this modified Van de Steene & Zijlstra sample 
(Set 18 in Table 1) is negligible. This again supports the 
idea that these sample sets are not strongly affected by 
the Malmquist bias. 



According to criteria defined in this paper, the corre- 
sponding slopes for the Bulge samples represent the result 
of scattering in the L — D plane. This slope ((3 ps 2) was ob- 
tained for the extragalactic samples of SNRs (except M82 
sample) where Malmquist bias also does not have severe 
effect. Again, this is because all of the SNRs in each set 
are at approximately the same distance (see lUrosevic 2 002 
and lUrosevic et al. 2005j) . 



By examination of sample sets 5 - 9 in Table 1 (and 
Fig. 3) we try to demonstrate the effect of Malmquist bias. 
For set 5 with distances less than 0.4 kpc, the correspond- 
ing £ — D slope \s (3 — 1.47. Further distance comparisons 
(in the form of distance interval in kpc followed by the 
E - D slope) include: 0.4 < d < 0.6, 2.18; 0.6 < d < 1.0, 
2.37; 1.0 < d < 2.0, 2.64; and > 2.0, 4.44. This implies 
that an increase in distance gives a £ — D slope that is 
also greater. 

Using both the Tukey-Kramer (T-K) and the Schcffc 
(S) methods, we found the £ — D slope for sample set 
9 in Table 1 (> 2.0, 4.44) is statistically different at a 
confidence level of a = 0.05, if compared to the slopes 
of sets 5, 6, and 7. However, since all of these samples 
consist of PNe with some uncertainty in their distances, 
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any firmer conclusion about their Malmquist bias should 
be avoided. Using the same methods, sets 1 (< 1, 2.61) 
and 3 (> 2, 3.96) in Table 1 also represent a statistically 
different population at a confidence level of a = 0.05. The 
difference between their corresponding £ — D slopes is 
statistically significant at a confidence level of a = 0.2 
and can be explained by the Malmquist bias. Sample set 
2 in Table 1 is omitted from our slope analysis because of 
its extremely high level of scattering in the L — D plane 
(see Table 1, r L - D = 0.08). 

Inspection of Table 2 reveals that while all £ — D slopes 
are relatively close to the trivial form, set 4, 6 and 7 are a 
bit steeper. The increased slope of set 4 can be explained 
by perturbations induced by twelve very luminous and 
highly £ — D plane dispersed "expansion" PNe added to 
sample set 2. Since the gravitational and spectroscopic 
methods are mostly applied for distant PNe, slopes in sets 
6 and 7 may be steeper due to a Malmquist bias. However, 
all of the L — D correlation coefficients in Table 2 are very 
low and therefore poor for statistical investigations. 

The high level of scattering in the L — D plane, as seen 
in Figs. 1 (right) and 4, supports the idea that most of the 
slopes in Tables 1 and 2 do not have any real physical inter- 
pretation. It is this luminosity-diameter scattering artifact 
which produces the trivial £ oc D~ 2 form. Therefore, 24 
of these trivial £ — D slopes are not useful for the deter- 
mination of distances. Causes for this scattering include: 
imprecisely determined calibrator distances, mixtures of 
different PN types in the same sample, limitations in sen- 
sitivity and resolution of radio surveys, source confusion 
and to a lesser extent, Malmquist bias. 

Although, our criteria imply that the relation given in 
Equation (30) could be useful for distance determination, 
we assign a low confidence that it represents a real evolu- 
tionary track. The £ — D slope from this relation repre- 
sents the result of coupling of evolutionary and selection 
effects. The corresponding sample set consist of only 13 
PNe and they do not belong to the same morphological 
type. Thus, this sample is not complete or homogeneous 
and distance estimates obtained by using Equation (30) 
should be viewed with appropriate caution. 

PN radio surface brightness depends on a number of 
factors that include: intrinsic nebula evolution, changes 
in density and temperature, changes in medium ioniza- 
tion and extremely large luminosity variations in the cen- 
tral star during its evolution. In Section 2, we assumed 
that shell expansion was constant throughout each phase 
of nebula evolution. On the contrary, PN shells are ac- 
celerated by the central star fast winds at the beginning 
of PN evolution, and only later are they driven at ap- 
proximately constant speed by the hot bubbles created by 
the crossing of the fast winds through the inner shocks. 
Shaw et al. (2001) discusses acceleration of PN shells, i.e. 
how the velocity of PN shells is related to the evolution 
of the central star. They emphasized that the accelera- 
tion of shells is nearly zero in larger (and probably older) 
PNe having diameters > 0.03 pc. Almost all of the PNe 
from sample sets presented in Tables 1 and 2 are larger 



SSV (d<1) 
SSV(1 <d<2) 

SSV (d > 2) 
d < 1 kpc 
d > 2 kpc 




D[pc] 

Fig. 3. The £„ - D diagram at 5 GHz for all 39 SSV 
PNe (sample set 4) with reliable individually determined 
distances. The straight lines represent the least square fit 
lines through the two distance limiting Galactic samples. 
The difference in slope of corresponding £ — D relations 
are defined for PNe with distances < 1 kpc and > 2 kpc 
(represented with dash and tick line, respectively). These 
could be explained as a result of the volume selection effect 
called Malmquist bias. 



than 0.03 pc (see Fig. 4) and therefore in later phases of 
evolution where acceleration (or central star influence) is 
negligible. This explanation is admittedly simplistic. For 
essential understanding of the theoretical and empirical 
£ — D dependencies, the full explanation of the complex 
co-evolution of a nebula and its central star would be 
needed. 
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The most massive central stars will never fully ionize 
their environment and those PNe will never become op- 
tically thin at 5 GHz since their central stars evolve so 
rapidly. Those observed PNe with smaller diameters and 
surface brightnesses (optically thick PNe at 5 GHz) are 
located significantly lower in the £ — D or L — D planes 
when compared to the average (thin) PNe. Their evolu- 
tion are surely not defined by the £ — D relations derived 
here. Although the majority of the thick PNe will become 
thin, the question is WHEN? This depends on a number of 
factors during their evolution. These effects significantly 
contribute to uncertainties in addition to the huge scat- 
tering found in the L — D plane. 

Other factors not included in our analysis may have 
a large impact on the evolution of PN emissivity. These 
include: compression of gas induced by the interaction of 
fast wind with the medium on the main shell inner bound- 
ary, isothermal shocks associated with the passage of D- 
ionization fronts through neutral gas and progressive ion- 
ization of greater amounts of surrounding neutral gas (e.g. 
IShaw et al. 20061 and references therein). We do not an- 
alyze ionization fronts in this paper; our derivations in 
Section 2 are based only on the dynamics of correspond- 
ing shock waves. The surface brightness evolution may be 
more closely connected to ionization than the dynamical 
expansion of the shell. This may be good enough reason to 
challenge the quality of our derived theoretical relations 
and an additional explanation for scattering in the L — D 
plane. Finally, we concede that some PNe could have dif- 
ferent initial conditions leading to independent evolution- 
ary paths. These paths could follow the same theoretical 
£ — D curve but with varying intercepts, therefore leading 
to the observed scatter. 

In order to define a valid PN sample for statistical 
investigation, one needs to separate a higher number of 
nearby PNe with the same morphological characteristics. 
We define such a sample (USNO-PN) in the subsection be- 
low (Section 5.1). However, the quality of radio-continuum 
observations today simply is not sensitive enough to pro- 
vide a high number of similar PNe having reliably deter- 
mined distances. 



5.1. The USNO-PN Sample 

We have examined the morphological properties of 
PNe from the USNO-PN sample using the radio im- 
ages of these objects found in the literature. We use 
comparisons with radio morphological classes given by 



Aaquist & Kwok (1996) which are based on the prolate 



ellipsoidal model (PES). The PES model describes diver- 
gency in observed (empirical) morphologies through dif- 
ferent sky projections of relatively simple spherical shell 
models with both radial and latitude density gradients and 
with different ionization depths. USNO-PN sample objects 
can be classified as follows: circular (NGC 7293 or Helix 
Nebula, see |Zijlstra et al. 1989D ; open elliptical (NGC 
6853 or Dumbbell Nebula, see IBignell 19831); elliptical 



(NGC 6720 or Ring Nebul a, see|George et al. 1974 P ; sym- 
metric (A21 or S274, see lSalter et al. 1984 ft ; and S-tvpe 
(A243). However, a further search of the literature reveals 
that the intrinsic morphologies for all five of these objects 
is most likely to be bipolar, i.e. that the structures of the 
emitting regions (or shells) range from a thick disk (Helix 
Nebula) to a triaxial ellipsoid (A24, and Ring Nebula) 
to a barrel-like shape (A21, and Dumbbell Nebula) 
(see IO'Dell"T998l [Hen Fy et al. 19991 IHua fc Kwok 19991 
IQ'Dell et al. 20071 and IMeaburn et al. 2000p . The sixth 
object for which we found reliable radio observations (A7) 
appears to be a "classical" PN with a well defined spher- 
ical shell (jXilouris et al. 1996jK 

and 



We 



note 



that 



Zuckerman fc Gatley (1988) 



Kastner et al. (1996) 



detected a shock-excited 
from the Helix, Ring and Dumbbell 
an ionization bounded case for 
spectroscopic evi- 



H2 emission 
nebulae, implying 
these objects. Additionally, no 
dence for the existence of a high-velocity stellar wind 
has been found for the previously mentioned ob- 
jects, including A7 (|Cerruti-Sola fc Perinotto 1 985. 
IPatriarchi fc Perinotto 1991|) . Central stars of these ob- 
jects are most likely well pass the hydrogen-shell burning 
phase and approaching the white dwarf cooling sequence. 

A wide range of shell structures is predicted in the 
hydrodynamical models based on the ISW model and it 
is very likely that a large number of PNe classified by 
apparent morphology as spherical or elliptical, may in 
fact posses a bipolar structure. Even though morphologi- 
cal classifications based on optical imaging favor elliptical 



objects, Zuckerman & Aller (1986) predict that approxi 



matcly 50% of all PN in the Galaxy arc actually bipolar. 
Nevertheless, we conclude that the USNO-PN sample is 
not representative of the majority of PNe because of its 
"narrowness" in morphology (i.e. almost all appear to be 
bipolar). 

6. Summary 

The main results of this paper may be summarized as 
follows: 

(i) We have derived a new theoretical £ — D relation 
for PNe in the form £ oc D^ 1 by including the in- 
teraction between the AGB star wind and the fast 
wind from the central star. This dependence is ob- 
tained for both momentum conserving and energy 
conserving phases of evolution. Only for the final 
stage, we did derive the same theoretical dependence 
(£ cx D~ 3 ) as in Paper I. 

(ii) We derived empirical £ — D relations for 6 updated 
PN samples from the literature, 12 new sample sets 
extracted from these and 9 additional sets as defined 
in this paper. We discuss the selection effects that in- 
fluence these PN sample sets. Our results show that 



7 For A24 we acquire a 1.4 GHz radio image from the 
NRAO VLA Sky Survey (NVSS) postage stamp server 
(http: / /www.cv. nrao.edu/nvss/postage.shtml ) . 
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updated Galactic PN samples do not severely suf- 
fer from volume selection effects, i.e. Malmquist bias 
(similar as in case of the extragalactic SNR sam- 
ples). We do note that the E — D slope does increase 
with sample distance, possibly due to Malmquist 
bias as shown on two statistically different SSV sam- 
ples with distance limits < 1 and > 2 kpc. Because 
of this trend, we feel that samples consisting only of 
nearby PNe are useful for future statistical analysis, 
(iii) From the analysis of L — D dependencies presented 
here, we conclude that 24 empirical S — D relations 
(listed in Tables 1 and 2) for Galactic PNe should 
not be used to determine distances to other observed 
PNe or to establish any other PN evolutionary paths. 
These relations represent only a general trend of 
changes in surface brightness. Two obtained steep 
£ — D slopes for distant samples (d > 2 kpc) are 
probably the result of Malmquist bias. However, we 
extracted one useful nearby SSV subsample (d < 1 
kpc, set 1 in Table 1). This updated sample consists 
of PNe with reliable and individually determined dis- 
tances. A reasonable L — D dependence does exist for 
this nearby sample. Since this sample consists of only 
13 morphologically different PNe, the corresponding 
£ — D relation for distance determination must be 
used with caution. 
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Fig. 4. The L v — D diagram for the samples of PNe given in Table 1. 



